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SYNOPSIS 


Ultrasonic means have been widely used in a variety of engineering applications for the 
nondestructive evaluation of pipes and tubes. Defects such as cracks and inclusions act 
as a source locations of wave scattering when illuminated by an incident ultrasonic pulse 
through reflection, diffraction and mode-conversion. Interaction of elastic wave with cracks 
provide all the necessary information regarding the inverse characterization of the defects 
which has not yet been thoroughly solved yet. The study of forward scattering problem 
or the determination of wave field scattered by defects would provide more knowledge in 
implementing the inverse algorithms. Analytical solution for the forward scattering is avail- 
able only for particular scattering solutions. 

There are several popular NDE methods that have their own advantages and disadvan- 
tages for the practical inspection work. Recently, Ultrasonic guided waves have been found 
to have some specific advantage in tubing inspection. The doininent benifit is the pene- 
tration ability to inspect the entire length of the tube so that the ultrasonic guided wave 
technique can quickly locate the defects in tubing without moving the ultrasonic transducer. 

In the present work a finite element model has been developed to study the elastic wave 
propagation and scatering in thin tubes. Interaction of the guided wave with the cracks 
has been studied, with the basic aim of detecting the axial as well as the angular location 
of the crack. 

Prom the results it is concluded that a combination of three source locations circum- 
ferentially located at an angle of 120 degrees aparts is very effective in scanning the entire 
surface of the tube. The loading in the tangential plane proved effective. Some guide lines 
are provided to locate the axial and circumferential position of the cracks. 


© © © 
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Chapter 1 
Introduction 


1.1 Introduction 

Understanding the propagation of waves through elastic media is of potential benefit 
in several fields of science and engineering. The effects of very short term impact loads on 
the structures can be understood fully if analyzed in terms of propagation of stress waves. 
Crack propagation or the interaction of dynamic stress waves with the existing cracks, voids 
or inclusions in a material can be predicted if the dynamic effects are taken into account. 
Further the knowledge of elastic wave propagation is essential in a host of other fields such 
as ultrasonics, geophysics, electronics etc. 

Ultrasonics is a major field of application of wave phenomena. The interest here lies 
in the ultrasonic nondestructive evaluation of materials which is widely used in aerospace, 
transportation, chemical and energy industries for the detection and characterization of 
defects such as cracks and inclusions in critical structural members. Such tests are needed 
both in the manufacturing phases of the products in order to assure the quality as well as 
during their operating life time in order to prevent failure and asses the residual life. 

In essence, ultrasonic guided wave inspection involves introducing a stress pulse into 
the material under inspection and observing the defects such as cracks and inclusions act 
as source locations of wave scattering when illuminated by an incident pulse through reflec- 
tion, diffraction and mode-conversion. An adequate understanding of the wave propagation 
is essential for the optimization of the current NDT techniques and also for the develop- 
ment of the new techniques. 
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1.2 Literature Survey 

Elastic wave propagation and scattering has been under investigation in many fields 
such as seismology, geophysics, electronics and applied mechanics for many years. An ex- 
cellent review of the state of art of the theory of stress wave propagation and scattering 
was given by Pao[l]. 

The current developments in high speed computers, range of work stations and graph- 
ics packages are all bringing revolutionary changes in the numerical modeling of ultrasonic 
NDT. Examples'of most modeling techniques employed to study ultrasonic NDT have been 
reported in the series ’’Review of progresses in quantitative NDE” [2] and various other 
publications^, 4, 5]. 

1.2.1 Guided waves and NDT 

There are several popular NDE methods, such as bulk wave ultrasonic, radiography acoustic 
emission, eddy currents and visual-optics, that have their own advantages and disadvan- 
tages for practical inspection work. Recently developed ultrasonic guided wave techniques 
have some benefits in tubing inspection. 

Guided waves, are nothing but waves propogating through bodies with narrow 
cross sections. Plate, or Lamb waves which propagate in thin layers of materials are some 
examples. These waves, unlike bulk waves, excite the whole cross section of the thin struc- 
ture and propagate in a specific direction. This specific advantage of the guided waves has 
lead to its use in NDT for defect classification and characterization. 

The dominant benefit is the penetration ability to inspect an entire length of the tube so 
that the ultrasonic guided wave technique can quickly locate the defects in tubing without 
moving the ultrasonic guided wave transducer. 

Guided wave propagation along the pipe provide a potentially very attractive inspection 
technique. They may be excited from the pipe outer wall and say, only a key hole in the 
insulation can be made to allow access for a transducer system. The waves can propagate 
tens of meters along the pipe, so a long length pipe may be inspected from each transducer 
position. 

Many workers have investigated the use of guided waves for the inspection of plate like 



structures and there has been a considerable amount of work on there use in the inspection 
of pipes and tubing[23-26]. 

The first general solution of harmonic wave propagation in an infinitely long elastic 
hollow cylinder was obtained by Gazis and his coworkers using elasticity theory (27-28]. 
Armenaks, Gazis and Herrmann suminerized their work into a book entitled ’’Free Vibra- 
tions of Circular Cylindrical Shells” including numerical results of some ultrasonic guided 
wave modes [29]. Ultrasonic guided wave propagation in hollow cylinders was studied by 
Fitch[30]. For the generation of ultrasonic guided waves in heat exchanger tubing, Silk 
and Bainton[31] investigated piezoelectric ultrasonic probes to access the inside layer of the 
tube. Alleyne et al. [32] suggested the use of the fastest mode in a non-dispersive region 
of frequency range to minimize dispersion effects over long distances. For defect detection 
and optimization, Dirti et al. [33,34] discussed mode selection criterion based on the acous- 
tic field across the thickness of the of the wave guide, explaining why certain modes are 
more sensitive to certain geometries of defects than others. Rose et al. [35,30,37] employed 
the mode control concept and reported experimental results for multiple crack detection, a 
water insensitive mode, and long range ^inspection. There is some research in defect char- 
acterization and sizing studies using ultrasonic guided wave technique in both experiment 
and theory [38, 39]. Very little has been done for the analytical solution of the acoustic wave 
field in tubing for an applied non-axisymmetric loading. Dirti arid Rose[40] published a 
paper which contained the normal mode expansion technique to obtain the acoustic field 
distribution along a tube for an arbitrary applied surface traction. 

Alleyne et.al have worked on a project whose ultimate aim is to develop a guided wave 
testing technique for the inspection of pipe work in chemical plant, the target being to 
detect any areas of corrosion larger than 3Tx3T in area and T/2 deep where T is the 
pipe wall thickness. The technique is to work on a insulated pipe in the 2-12 inch nominal 
bore range and an inspection range of at least 15m from the transducer position is required. 

A first concern on the use of ultrasonic guided waves in tubing is to find the useful 
modes. Theory of elasticity and basic wave mechanics can be used to predict the infinite 
number of modes that can be generated in a tube. The inodes can be represented on 
dispersion curves for given frequency and phase velocity ranges. Also the group velocity 
dispersion curves can be used for an identification of received modes, mode selection, and 
defect location. 
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Guided wave field analysis or wave structure through the thickness of tubing is one of 
the impOl ln.nl, piMninci.cn: in selecting inodcd lot a i n lain pmpunc. It in shown in a piLpci 
by Pilarski, Ditri, and Rosc[25] that points of dominant in plane displacement are achieved 
for a phase velocity close to the longitudinal wave velocity. Points in the neighbourhood 
of these special points have dominant in plane components but still some out of plane 
component which allows coupling and minimal leakage of energy to fluid loaded tubes from 
the outside if this situation arc encountered. From a defect detection point of view, the 
displacement field analysis gives a good explanation of the sensitivity changes to defects. 
For long range screening purposes, it is better to use inodes with less outside sensitivity 
or/and lower frequencies. Modes with power concentration on the surface are more sensi- 
tive to surface cracks than the modes with power concentration on the center line of the 
tube thickness. 

Great progress is being made on the use of axisymmetric ultrasonic guided waves in 
tubing for defect detection, classification and sizing. Benefits of excellent penetration and 
sensitivity to certain defects are possible. However, ultrasonic guided wave propagation and 
reflection in tubing is still not completely understood. For more comprehensive knowledge 
of ultrasonic guided wave tubing, a study on non-axisymmetric modes cannot be avoided. 
A number of non-axisymmetric modes close to the axisymmetric modes on a dispersion 
diagram often makes it difficult to generate a single axisymmetric mode in tubing. Also, 
non-axisymmetric reflections occur due to non-nonaxisymmetric defects, even for axisym- 
metric impingment. 


1.2.2 Use of Finite Element Method in Wave Propogation Anal- 
ysis 

The development of appropriate method.*: for the aualysiH of shell structures is in 
creasingly demanded to ensure the integrity of structural designs. Analytical solutions to 
shell structures are limited in scope and in general are not applicable to arbitrary shapes, 
load conditions, irregular stiffening and support conditions, cutouts and many other as- 
pects of practical design. The finite element method has consequently become prominent 
in the analysis of such shells in view of the ease with which such complexities can be dealt 
with. 

Although finite element procedures have been applied to shell analysis for over twenty 
years, the search for suitable elements which are readily applicable to general shell struc- 
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tures is still in progress[6-9]. Tlireo distinct approaches to the finite element representation 
of shell structures have been employed based on 

(i) flat triangular or quadrilateral elements, 

(ii) curved elements formulated on the basis of various shell theories and 

(iii) elements derived from three-dimensional elements by the use of the degeneration 
methods. 

The first numerical studies of shell problems involved facetting the shell by plane trian- 
gular elements, onto which a membrane stiffness was superimposed . The results obtained 
were found to be satisfactory, but depending on the problem, very fine meshes had to be 
employed. A number of difficulties and shortcomings arise from the application of flat 
elements to curved shellsflO], such as the presence of ‘discontinuity’ in bending moments, 
which do not appear in the continuously curved actual structmes. 

In the second approach a classical concept is employed, whereby a shell theory is used 
as the starting point of the finite element formulation. Different shell theoi'ies have been 
developed from the three dimensional field of equations by incorporating various assump- 
tions appropriate to the structural behaviour. A variety of finite elements with different 
degrees of complexity have been formulated for deep and shallow shells. Curved shell el- 
ements based on the Kirchhoffs-Love theory guarantee a high solution accuracy but are 
complicated by convergence and compatibility requirements. In the case of a fully nonlin- 
ear analysis, additional difficulties arise in the formulation of efficient finite elements due 
to the unavalibility of a general nonlinear shell theory. 

The third approach employed in the formulation of shell elements avoids the complexi- 
ties of fully general shell theories by discretizing directly the three-dimensional equations of 
continuum mechanics. Isoparametric elements with independent rotational and displace- 
ment degrees of freedoin(dof) are employed, in which the three dimensional stress and 
strain conditions are degenerated to shell behaviour. The procedure was originally intro- 
duced by Ahmad et al. [11] for the linear analysis of moderately thick shells. Although 
this type of element seemed to be very promising when it was introduced, difficulties later 
appeared, due to the degeneration process, as the thickness of the element was reduced. A 
great improvement of the model was achieved by means of the so-called reduced integration 
technique. 
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In numerical modeling applied to elastic waves much effort has been concentrated on 
the finite difference method and more recently, on the finite element method. These models 
consider mid-frequency scattering where complex mode-conversion phenomenon is encoun- 
tered. Such models give a complete description of the time development of the wave-field 
of interest. They can also model mixed solid-fluid problems and can model a wide range of 
material properties. 

Developments in the finite difference modeling of ultrasonic wave propagation were re- 
ported by Bond et al [15]. Results were obtained for a range of two dimensional geometries 
with pulsed plane wave excitation and found to agree with experiments. In addition immer- 
sion inspection systems were modeled where there is propagation in both an elastic solid 
and a viscous fluid. Finite difference model results were produced for pulsed wave scatter- 
ing by a semi-infinite 2-D slit and compared with those obtained from GTD[16]. Bond & 
Saffari[17] has emphasized the mode-conversion technique in ultrasonic NDT and has given 
experimental results of crack detection by making use of the mode-conversion phenomenon. 

The use of FEM to model ultrasonic scattering is rather resent. Its superiority lies in its 
ability to handle complex shaped domains and defect shapes, discontinuities, inhomogen- 
ilies, nonlincarities in material etc. The pioneering effort in this direction has been by 
Lord [20,21] and his co-workers who produced a scheme which employ a finite element for- 
mulation for spatial discretization and finite difference scheme for temporal discretization. 
Ludwig and Lord [20] reported details of a 2-d finite element code employing Newmark’s 
scheme for time integration. Results in the form of A-scans and visualization of wave fields 
were produced for forward scattering by a circular hole in a homogeneous isotropic medium. 
Transmitter and receiver modals were also incorporated. You and Lord [21] reported im- 
provements from computational point of view, by employing a central-difference formulae 
for time integration and using a lumped mass matrix. They also studied wave interactions 
with a rectangular slot and surface breaking cracks in isotropic solids. 

A detailed analysis of the various time integration schemes from the point of view of 
accuracy and speed, was reported by Seron et al [22]. It was concluded that, even though 
central-difference scheme was the best in terms of computational speed and cost, New- 
mark’s constant average acceleration scheme gives better accuracy. 
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1.3 Present work and Thesis Layout 

In the present work a finite element model to study elastic wave propagation and 
scattering in thin tubes with cracks has been developed. A methodology has been suggested 
to for the location and classification of the defects. 

Chapter 11 of the thesis gives the basics of the guided wave propagation in thin layers. 


Chapter 111 describes the finite element modeling of the wave propogation in shells. The 
problem of spatial and temporal discretization have been discussed. 

Chapter IV presents the results for the validation of the computer code and those of 
the simulations of the tubes with cracks. . The conclusions and the scope for future work 
are given in chapter V. 



Chapter 2 


Basics of Guided Wave Propagation 
In Thin Layers 

2.1 Introduction 

This chapter deals with a brief discussion of elastic wave propagation in the isotropic 
thin layers. For plates, exact solution of the equation of motion can be derived from 
the classical elasticity theory involving infinitesimal deformations. These solutions satisfy 
the equation of motion and the boundary conditions related to traction-free surfaces. In 
addition they represent the modes of propagation which exist in simple types of elastic 
wave guides. 

2.2 Propagation in Thin Layers 

2.2.1 Integration of equation of elastic motion 

The vector differential equation governing the small elastic motions of an isotropic 
elastic medium in the absence of body forces can be written as 

eft U 

pA 2 u + (A + /z)V(A.n) = P-Qfi (2-1) 


where 

^displacement 

t=time 

A,p=Lame’s constants 
p=density 
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In the particular problem of analyzing motions in the thin layers, solutions of Eq.(2.1) 
are to be found consistent with certain boundary conditions(b.c’s) to be proscribed. The 
analysis will lead to a discrete but infinite set of solutions for u which are associated with 
inodes of propagation. The solution of Eq.(2.1) can be recast in a combination of a vector 
potential function T and a scalar potential d> such that 

u = V<I> + (Vx£) (2.2) 


If $ and d/ are solutions to the equations 


V%V 2 < E> = 


a 2 $ 

dt 2 


(2.3) 




dH 

~d¥ 


(2.4) 


and 


V.£ = 0 


(2.5) 


where 


Vi = ( A - + 2 ^ )V 2 


V a = (-) 1/2 (2.6) 

P 

Solutions to the component equations will be sought corresponding to wave motions prop- 
agating in the positive axial direction. By assuming a set of solutions of the form 

<3> = (Acos(cta:) + B sin(a:a;))e t(7 * _ ‘ Ji) 


{(yz'-ut) 


= ( C cos (/3x) + Dsin(/3x))e‘ 
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= {E cos(/3x) + Fsin(Px))e i{ ' rz - wL) 


V, = (G cos(#e) + H sin(/3x))e i(72 “^ (2.7) 

Direct substitution of these potential function into Eq.(2.3) and (2.4) indicates that they 
are the solutions if 


a 2 + 7 2 = u 2 /V 2 

(2.8) 

P 2 + 1 2 =m 2 /V 2 

(2.9) 


Substituting these equations in the displacements obtained from Eq(2.2) and then us- 
ing the b.c’s and the stress strain relationships six equations in A,B,C,D,E,F,G, and H are 
obtained. The other equations from Eq(2.5). 

These eight equations can be considered as a set of homogeneous linear equations in 
the eight variables A-Ii. A necessary and sufficient condition for these equations is that the 
determinant of the coefficients of the variables be zero. 

The determinant formed by the these set of equations can be written as a product 
of four subdeterminants. The equation have solutions for four separate and independent 
conditions; namely, when any one of the four subdeterminants is equal to zero. 

Solution I. 

D,G^0 

XL's — 0 

u y =((3G + iD){ sm/3x)e^ z -^ (2.10) 

u z = 0. 


Solution II 
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C,H^() 

u x =0 

u v = i~/3H + iC)(c os px)e i( ' z ~ wt ) 
u z = 0. 

Solution III. 

A, F^0 

u x = (— olA sin ax — iF sin fix)e i( - z ~ wt ' ) 
u y =0 

u z = (—fiF cos fix — iF cos ax)e l ^ z ~ wt ^ 

Solution IV. 

B, E#0 

u x = ( aB cos ax — iE cos fix)e l ^ z ~ ult ) 
u y =0 

u z — (—fiE sin fix + iB sin ax)e^ z ~ wt ^ 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


For each of the four solutions, the equation formed by setting the appropriate subde- 
terminant equal to zero is called a frequency equation. 

In the discussion of wave motion in plate, it will be convenient to classify the solutions in 
terms four families of wave motions. In the case of both Solutions I and II, the displacement 
is transverse to the propogation direction of the wave. Thus, Solution I and II corresponds 
to shear waves with a particle displacement parallel to the boundary surfaces. Waves of this 
type are often called SH waves and may be described in terms of antisymmetric SH waves 
with particle displacements of the form given by Eq(2.9) and symmetric SH waves the form 
given by Eq(2.10). Fig(2.1) gives the frequency spectrum of symmetric and antisymmetric 
shear modes in an infinite plate. It is observed from the Fig(2.1) that for all symmetric 
and anti-symmetric modes the group velocity and phase velocity are functions of frequency, 
indicating that these velocities become imaginary for values ujb/V 3 $. Free propogation in a 
given mode does not occur at frequencies lower than ub/V s = r<3>. At this point the group 
velocity in a given mode becomes zero and the frequency at which this happens is called 
the cut-off frequency of the mode. 

The wave motions associated with Solutions III and IV are somewhat more complicated 
than those associated with the SH modes: first, in the sense that the two displacements 
components are nonzero, and second, the total wave motion involves a combination of shear 
and dilatation waves . In Solution III, the particle displacement vectors are symmetrically 
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arranged with respect to the median plane; so that this is asset iated with the symmetrical 
or longitudinal family of modes. The particle displacement vectors in Solution IV are an- 
tisymmetrically arranged with respect to the median plane arid this solution is associated 
with the antisymmetrical or flexural family of modes. Fig(2.2) gives the frequency spec- 
trum of longitudinal and flexural modes in an infinite plate. 


2.3 Guided waves in tube inspection 

There are several popular NDE methods, such as bulk wave ultrasonic, radiography acoustic 
emission, eddy currents and visual-optics, that have their own advantages and disadvan- 
tages for practical inspection work. Recently developed ultrasonic guided wave techniques 
have some benelils in tubing inspection. The dominant benefit is the penetration ability 
to inspect an entire length of the tube so that the ultrasonic guided wave technique can 
locate the defects at any point in tubing without moving the ultrasonic wave transducer. 
The concept of using guided waves in tubing inspection has been published for many years. 

A first concern on the use of ultrasonic guided waves in tubing is to find useful inodes. 
There are at any frequency, at least two guided wave modes available and the number in- 
creasess as does the inspection frequency or layer thickness. Theory of elasticity and basic 
wave mechanics can be used to predict the infinite number of modes that can be generated 
in a tube. It is known that, two families of Lamb waves exists; the symmetric and the 
antisymmetric, and the relationship between phase velocity and product of frequency and 
thickness for different families is given by, 


f2 s (w^, Vp/i) 

Symmetric Modes 


tan(fc ts 6/2) 4 (Pkgkts 
tan(fc a 6/2) {k? s - p 2 ) 2 


kp/i) 

Antisymmetric Modes 


tan(fej S &/2) (fcg, - /? 2 ) 2 _ 
tan {ktib/2) 4p 2 kuk ts 


(2.14) 


(2.15) 


where k u and k ts are given by, 
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4 = (ff) 2 - P 2 

4, = (ft) 2 - P 2 

The wavenumber, /?, is numerically equal to w/V p } x where V p h is the phase velocity of the 
Lamb wave mode and to is the circular frequency. The phase velocity of a wave is related 
to its wavelength by the relation,!^ = ( co/2tt)\ . The quantities Vi and V s , which represent 
the longitudinal and shear wave velocities in the bulk material, contains the information on 
the material properties oi the material in which the Lamb waves propagate. The velocities 
are given by, 

Vi = ( A ^ 2 / V /2 (2.16) 

and 

V, = (~) 1/2 (2.17) 

P 


where A and p represent the Lame constants of the layer and p represents its mass 
density. 

For a given frequency, Eqn(2.14) and Eqn(2.15) can be considered implicit transcedental 
functions of phase velocity alone, and they will be satisfied by an infinite number of real, 
imaginary or complex values of phase velocity. The real solutions for a given frequency 
represent undamped propogating modes of the structure whereas the imaginary and com- 
plex roots represent exponentially decaying modes which do not propogate. 

The symmetric Lamb wave modes, governed by Eq.(2.13), have deformation-fields which 
are symmetric about the midplane of the layer, while the anti-symmetric modes, governed 
by Eq.(2.14) have deformation fields which are anti-symmetric with respect to the midplane 
of the layer. 

For given parameters A, p and p and a given frequency u times plate thickness (d) 
product, the real roots (/?) of Eqs. (2.13) and (2.14) represent undamped propagating 
waves permissible in the layer. There are, at any given wd a finite number of such roots, 
with a minimum of two at n>d =0 and increasing as to increases. Snell’s Law provides the 
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relationship between the angle of incidence and phase velocity of the guided waves. 


Qi = sin 


V T 


ph 


(2. IS) 


where V[ represents the velocity of longitudinal waves in the medium containing the 
incident waves. In practice, where finite sized transducers are used, there is actually a 
phase velocity spectrum which is excited by the incident wave, not the single phase ve- 
locity implied by Eq.(2.17). The range of phase velocities which may be excited by a 
given transducer depends upon the size of the transducer relative to the wavelength of the 
generated inodes as well as on the pressure distribution across I, he face of the transducer. 
Fig(2.3) shows the phase velocity dispersion curve for aluminium. This dispersion curve 
gives information about the modes that propogate in the plate at the particular frequency, 
thickness product. The larger the transducer, the more selective it is to the ’’Snel’s law” 
phase velocity given by Eq.(2.17). 

Lamb waves propagate at their group velocities, V g , as opposed to their phase velocity, 
Vph, and the two velocities are not, in general, equal. The group velocity dispersion curves 
for the symmetric or antisymmetric Lamb wave modes can be calculated from the appro- 
priate phase velocity curves by the relation, 


r _ duj _ —d£l/d/3 
3 dp 00,/du) 


(2.19) 


■ where Q represents il s for the symmetric modes or H a for the antisymmetric modes. 

Guided wave field analysis or wave structure through the thickness of tubing is one of 
the important parameters in selecting modes for a certain purpose. It is by Pilarski et 
al[25], that points of dominant in plane displacement are achieved for a phase velocity close 
to the longitudinal wave velocity. Points in the neighbourhood of these special points have 
dominant in plane components but still some out of plane component which allows coupling 
and minimal leakage of energy to fluid loaded tubes from the outside if this situation are 
encountered. From a defect detection point of view, the displacement field analysis give a 
good explanation of the sensitivity changes to defects. For long range screening purposes, 
it is better to use modes with less outside sensitivity or/and lower frequencies. Modes with 
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FIG 2.1: FREQUENCY SPECTRUk OF SYMMETRIC AND ANTISYMMETRIC 
SHEAR MODES IN AN INFINITE PLATE. 
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FIG 2,2: FREQUENCY SPECTRUM OF LONGITUDINAL AND FLEXURAL MODES 
IN AN INFINITE PLATE. 
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Phase and group velocity dispersion curves for aluminum single layer. 


Chapter 3 

Finite Element Modeling of Shells 


3.1 Introduction 

The development of appropriate methods for the analysis of shell structures is in- 
creasingly demanded to ensure the integrity of structural designs. Analytical solutions to 
shell structures are limited in scope and in general are not applicable to arbitrary shapes, 
load conditions, irregular stiffening and support conditions, cutouts and many other as- 
pects of practical design. The finite element method has consequently become pro min ent, 
in the analysis of such shells in view of the ease with which such complexities can be dealt 
with. 

Although finite element procedures have been applied to shell analysis for over years, 
the search for suitable elements which are readily applicable to general shell structures is 
still in progress. Three distinct approaches to the finite element representation of shell 
structures have been employed based on 

(i) flat triangular or quadrilateral elements, 

(ii) curved elements formulated on the basis of various shell theories and 

(iii) elements derived from three-dimensional elements by the use of the degeneration meth- 
ods. 

The first numerical studies of shell problems involved facetting the shell by plane trian- 
gular elements, onto which a membrane stiffness was superimposed . The results obtained 
were found to be satisfactory, but depending on the problem, very fine meshes had to be 
employed. A number of difficulties and shortcomings arise from the application of flat ele- 
ments to curved shells, such as the presence of ’discontinuity’ in bending moments, which 
do not appear in the continuously curved actual structures. 
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In the second approach a classical concept is employed, whereby a shell theory is used 
as the starting point of the finite element formulation. Different shell theories have been 
developed from the three dimensional field of equations by incorporating various assump- 
tions assumptions appropriate to the structural behaviour. A variety of finite elements with 
different degrees of complexity have bcen fonmilated’for deep and shallow shells. Curved 
shell elements based on the KirchhofFs-Love theory guarantee a high solution accuracy but 
are complicated by convergence and compatibility requirements. In the case of a fully non- 
linear analysis, additional difficulties arise in the formulation of efficient finite elements due 
to the unavalibility of a general nonlinear shell theory. 

The third approach employed in the formulation of shell elements avoids the complexi- 
ties of fully general shell theories by discretizing directly the three-dimensional equations of 
continuum mechanics. Isoparametric elements with independent rotational and displace- 
ment degrees of freedom(dof) are employed, in which the three dimensional stress and strain 
conditions are degenerated to shell behaviour. Although this type of element seemed to be 
very promising when it was introduced, difficulties later appeared, due to the degeneration 
process, as the thickness of the element was reduced. A great improvement of the model 
was achieved by means of the so-called reduced integration technique. 


3.1.1 Element Formulation 

General 

Fig(3.1) shows a solid three-dimensional element based on a quadratic displacement 
field and and the corresponding quadratic degenerated shell element. Two basic assump- 
tions are adopted in this process: Firstly it is assumed that, even for thick shells, ’normals’ 
to the middle surface remain practically straight after deformation. Secondly, the strain 
energy corresponding to stresses perpendicular to the middle surface is disregarded, i.e. the 
stress component normal to the shell mid-surface is constrained to be zero in the constitu- 
tive equations. 

Five degrees of freedom are specified at each nodal point point, corresponding to its 
three displacements and the two rotations of the ’normal’ at the node. The definition of 
independent rotational and displacement degrees of freedom permits transverse shear defor- 
mation to be taken into account, since rotations are not tied to the slope of the mid-surface. 
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Coordinate Systems 

In the formulation of degenerate shell elements different coordinate systems have to 
be used. The four coordinate frames employed in the present formulation arc as described 
below. 

Global coordinate Set -Xi This is a cartesian coordinate system, freely chosen, in 
relation to which the geometry of the structure is defined in space. Nodal coordinates and 
displacements, as well as the global stiffness matrix and applied force vector are referred 
to this system. ■ 

The following notation is used Fig(3.1): 

Xi(i = 1, 3) and xi = x, x^ = y,x-s = 2 

Ui(i = 1, 3) and U\ = u, u 2 = v, — w 

Xi(i = 1, 3) is a unit vector in thea;; direction. 

Nodal Coordinate set -V ik 

A nodal coordinate system is defined at each nodal point with origin at the ref- 
erence surface (mid-surface). The vector V^ k is constructed from the nodal coordinates of 
the top and bottom surfaces at node k, 

Ya = I?'-3S' (3.1) 

where x k — [xk Vk *kY 

The vector V_ lk is perpendicular to V^ k and parallel to the global xz-plane, 

va=vs,va=o.o,^=-ia 


(3.2) 
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where the superscripts refer to the vector components in the global coordinate system. 
The vector V 2fc is perpendicular to the plane defined by V lk and V^ k , 


V 


■2k 


V* k x y 


it 


(3.3) 


The unit vectors in the directions of V 1 fc ,y 2A . and y Jjt are represented by y Ifcj V_ 2k andV^*. 
respectively . • 

The vector V defines the direction of the ’normal’ at node k, which is not necessarily 
perpendicular to the mid-surface at k. Vectors V_ lk and V_ 2k define the rotations ( and 
fiik respectively) of the corresponding ’normal’. The advantage of this definition for the 
vector]/^ (not necessarily perpendicular to the shell mid-surface) is that, as a consequence, 
there are no gaps or overlaps along element boundaries. 

Curvilinear Coordinate Set - £,77,£ 

In this system £ and rj are two curvilinear coordinates in the middle plane of the 
shell element and £ is a linear coordinate in the thickness direction. It is assumed that 
£, 77 and £ vary between -1 and +1 on the respective faces of the element. The relations 
between the curvilinear coordinates (£, 77, £) and the global coordinates (x,y,z) are given 
by expressions (-), which later define the element geometry. It should be noted that the £- 
direction is only approximately perpendicular to the shell middle surface, since £ is defined 
as a function of V^. 

Local coordinate Set -x\ 

This is a Cartesian coordinate system defined at the sampling points where(in) ? 
stresses and strains are to be calculated. Various directions are obtained by using the re- 
lations given below, 



dx 


dx 


oi 


dn 

y, = 
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OS 

dz 
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dz 


os J 


- dr) _ 


(3.4) 
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.2 — 2L3 ^ 


(3.5) 


(3.6) 


The local coordinate system varies along the thickness for any ’normal’ with this vari- 
ation depending on the the shell curvature and variable thickness. 


Element Geometry 

The global coordinates of pairs of joints on the top and bottom surface at each node 
are usually input to define the element geometry. In the isoparametric formulation the 
coordinates of a point within the element are obtained by applying the element shape func- 
tions to the nodal coordinates, 

^ = + (3.7) 

Jt=l Z k= l 


where i=l,3 refers to the three global directions, 
n - is the number of nodes per element, 

N k (£, jj)(k — 1 , n) are the element shape functions corresponding to the surface £=constant, 
h k - is the shell thickness at node k, i.e. the respective ’normal’ length, 

£, 77 , ( - are the curvilinear coordinates of the point under consideration. 



Displacement Field 

The element displacement field can be expressed as, 
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(3.8) 


The contribution to the global displacements from a given node k is: 
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(3.9) 


Calculation of Strains 

In order to easily deal with the shell assumption of zero normal stress in the z - 
direction a z > = 0, the strain components should be defined in terms of the local system of 
axes x\. The five significant strain components are, 
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(3.10) 


where u , v and w are the displacement components in the local system x\. 

Calculation of Stresses 

Taking into consideration the assumption of zero stress in the direction perpendicular 
to the shell mid-surface (ay = 0), the five stress components in the local system are: 
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(3.11) 


Layered model 

If the shell is built up from a series of different materials, such that the material 
properties(and stress) are discontinuous functions of £, an appropriate integration through 
the thickness has to be carried out. Here layered approach is employed, wherein a mid- 
point rule integration scheme is adopted for each layer. Each layer contains stress points 
on its mid-surface. The stress components of the layer are computed at these stress points 
and are assumed to be constant over the thickness of each layer, so that the actual stress 
distribution of the shell is modeled by a piecewise constant approximation. The specifica- 
tion of layer thickness in terms of curvilinear coordinate, C> permits the variation of the 
layer thickness as the shell thickness varies. 

The stress resultants are obtained by integrating the corresponding stress components 
with respect to the thickness coordinate: 

-normal forces, 

rh /2 n . 

N x = I a x dz = h/2 cr l x ^C 


rh/2 JL . 

Ny = / cjydz = h/2 ]T a l y AC (3.12) 

J-h/2 i =1 




a x zdz = -h 2 /AY, a iC A C 

«=i 


-bending moments 
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-shear forces 


hL-- 


rh/2 _ n 

V ~ I /,/•; 0 V Z ^ Z ~ ~d\ l j\ ^ vfX'AC 


Mi 


xy 


rh/2 n 

= / . o~ X yZdz = -h 2 /4j2<7xyiC&C 

J-h/2 i=1 


rh/2 n 

Qx ~ J_ h/2 Txzdz — hj 2 T xz^C 


rh/2 n 

Qy — Tyzdz = k /2 Ty Z A(~ 

J-h/2 i=1 


( 3 . 13 ) 


( 3 . 14 ) 


where n (i=l,n) is the number of layers. 

In the present code the strain matrix B_ is calculated at the mid-surface of each layer: 
the element stiffness matrix K_ e and the internal force f e are thus defined as follows: 


K‘ = J j*' B T DBd(dA 


( 3 . 15 ) 


l e = J J +1 B T aJd(dA 


( 3 . 16 ) 


where 

/ dA— f J+l drj£ (integration on layer midsurface) 

and J is the determinant of the Jacobin matrix. This process of integration in the thick- 
ness direction is computationally more expensive, but is more appropriate for thick and 
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variable thickness problems in which 


Uk> variation of the local .system of axes, as well its 


the variation of the Jacobian matrix through the shell thickness must be taken into account. 


Element Mass Matrix 

Due to computational convenience and simplicity lumping of mass matrices is con- 
sidered. Te essential requirement of mass preservation is given as: 


J2 Mii = f pdQ (3.17) 

i 

where Mu is the diagonal of the lumped mass matrix M. 

The approximation used here for mass lumping is Diagonal scaling procedure. This method 
satisfies the essential requirement given by the above equation. 

In this method 


Mu = aMu 

with a so adjusted that Eq(3.-) is satisfied. By following this method for our present anal- 
ysis the mass at the end nodes and and at the central nodes will be 1/36 and 8/36 times 
that of the total mass respectively. 


3.2 Time Discretization 

The system of linear, second order, ordinary differential equations with constant co- 
efficients has to be integrated in the time interval (0,T) by discretizing the time variables 
as t n = n A t, 0 < ji! < N where A t — T/N. 

3.2.1 Central Difference Scheme 

In the central difference scheme it is assumed that 


l dU 

dt 


1 /t+Atrj _ t-A trj\ 

2AV 1 



(3.18) 
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2 l U + t+At U) 


Consider the equation of motion at time t, i.e., 


tfPU 

+ K t U= t F 


(3.19) 


Substitution of Eq.(3.10) into (3.11) gives 

(-^*0 ‘ +A ‘ U = ‘ F ~ {K " ' u ~ X M) ! ' At£/ (3 - 20) 

from which one solve for t+At U, given the initial conditions. 

It should be noted that the solution of t+At U is based on the equilibrium conditions at 
time t. For this reason this integration procedure is called a(f explici^integration method. 
It is noted that such an integration scheme do not require a factorization of the stiffness 
matrix in the step-by-step solution. So if the mass matrix is diagonalized the assembly of 
the mass and stiffness matrix are not required. This is the important advantage of central 
difference scheme from computational point of view, even though it is inferior to the im- 
plicit Newmark’s scheme which uses the equilibrium conditions at time t, + At. 

With a diagonal mass matrix above equation can be simplified as 

t+M Ui = l Ri{—) (3.21) 

mu 


where 
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(3.22) 
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3.2.2 Newmark’s Method 

The Newmark’s algorithm under the assumption of constant-average acceleration is 
as follows, 


t+At dU 

dt 


l dU_ 

dt 


+ [(1 - S) 
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~W 


+ 6 


t +At d 2 U 
dt 2 ^ 


l+M u = L u + 
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A t + [(1/2 - a) 


l d 2 U 
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t+AtQ2jj 

dt 2 


]At 2 


(3.23) 


where a and <5 are parameters that can be determined to obtain an unconditionally sta- 
ble scheme known as the constant average acceleration method, in which 5 = 1/2 and a = 1. 

In addition to the above difference equations, for solution of displacements, velocities 
and acceleration, the equilibrium equations at time t + At are considered: 


t+Ato2rr 

M— — + K i+At U = t+At F 

dt 2 


(3.24) 


Using the above two equations we obtain the following equations for t+At U, 


[M + ~K\ ,+A ‘y = ,+M F + M{‘U + M~ + ^f~] 


(3.25) 


It is seen that this is an implicit scheme and requires the factorization of the stiffness 
matrix. This is reported to be more accurate than the central difference scheme. 
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FIG 3.1: ELEMENT GEOMETRY AND COORDINATE SYSTEMS. 



Chapter 4 

Results and Discussion 


This chapter discusses the results of guided wave propagation in a few problems 
using FEM as discussed in earlier chapters. The emphasis is on its application as a Non- 
Destrective technique(NDT) to detect cracks. 


4.1 Validation of the code 

In order to verify the code an example of cylindrical barrel vault was considered. 
The geometry and physical details are given in Fig(4.1). The barrel is supported on rigid 
diaphragms and is loaded by its own weight. In the present analysis 16(4 x 4) degenerated 
shell elements are used. Fig(4.2) shows some comparitive results, obtained by using the 
present code with that of exact results. The results are almost identical with the exact 
ones. 

Fig(4.2a) gives the vertical deflection of the central section of the vault, Fig(4.2b) gives 
horizontal displacement of the support. 


4.2 Wave Propogation in Tubes 

4.2.1 Details of Specimen geometry and loading 

The domain modelled is a thin tube of outer diameter 18.9mm, 1.22mm thick and 
350mm long. The specimen geometry is shown in Fig(4.3). Two different coordinate 
frames are used in the analysis viz., cartesian coordinate system and polar coordinatesys- 
tem. Fig(4.3) also contains the the two coordinate frames viz., cylindrical(r,0,y) and carte- 
sian(X,Y,Z) coordinate frames with respect to the specimen. 
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For convenience, the output of the Finite element code are expressed in radial, circuin- 
liToufinI mid axial component-H an hIiovvii in Fig('1.3). 

As seen in the Fig(4.3) ultrasonic pulse is applied at three locations viz., Si, 52 and 5 3 
which are circumferentially 120 degrees apart. The loads Fi,F 2 and F 3 are applied in the 
tangential plane at the loading points. Load at each of the three source locations is at an 

angle of 20 degrees to the axial plane. The loading as shown in the Fig(4.3) is specified 

in cartesian coordinate system. The transducer is modelled with a cosine pulse of one full 
cycle. Fig(4.4) contains the input at the three source locations in the cartesian coordinates. 

Eight-noded degenerated shell elements have been used. For time integration central 
difference and Newmark’s time integration scheme is used. 

The properties of the material are taken as follows: 

E=69.0 Gpa, G=26.4 Gpa, =0.33, p = 1800^/m 3 

The compressional, shear wave velocities in the material can be found as 

Ci — 7.2Qkm/sec , C s = 3.80 km/ sec, Cr — 3.10 km j sec 

The transducer pulse is modelled as 
f(t) = (1 — cos(2IIf)) 0 < t < 12.5 psec 

= 0 t > 12.5 psec 

A frequency of 0.1 MHz is selected, because for this product of frequency and thickness 
only two modes of guided waves exists as seen from the dispersion curves for aluminium in 
Fig(2.3). For this frequency the corresponding wave lengths are 

A / = 67.5 mm, X s = 31.0mm, Xr = 29.0mm 

Table (4.1) gives details of the five simulations performed. 

First simulations are performed on a defect free specimen, later two specimens with cir- 
cumferential cracks and two specimens with longitudinal crack at various axial and angular 
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First simulations are performed on a defect free specimen, later two specimens with cir- 
cumferential cracks and two specimens with longitudinal crack at various axial and angular 
positions are considered. 


Table(4.1) 


Cases 

No.of elein 

No.of nodes 

Axial location 
of crack (mm) 

Angular location 
of crack(inm) 

1. Defect free 

720 

2208 


“ 

2. Circumferential crack 

720 

2211 

175 

165-195 

3. Circumferential crack 

720 

2211 

233 

210-240 

4. Circumferential crack 

720 

2211 

163.3-18G.6 

180 

5. longitudinal crack 

720 

2211 

128.3-151.6 

120 


In all the simulations performed the load is applied at three locations which are circumfer- 
entially 120 degrees apart. The load at all the source locations was in the tangential plane. 

The finite element output includes radial(u), circumferential v), axial displacements(w) 
at the three source locations, at every time step. 

4.3 Propogation in defect-free specimen 

This simulation is performed to determine the movement of the wave propogation in 
the specimen and to verify the the speed of the wave in the axial and in the circumferential 
direction. Fig(4.3) gives the details of the geometry of the specimen and loading. 

Fig(4.5) and Fig(4.6) show the circumferential and axial dispalcement field in the do- 
main at different times after the application of the pulse. The displacement fields are 
plotted as isometric views to bring out the features. As the loading is in the tangential 
plane, displacement field is very small in the radial direction and hence those plots are not 
shown. The evolution and propogation of the deformations are clearly seen. In particular, 
from Fig (4. 5) and Fig(4.6) it can be seen that the axial displacements with the selected set 
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of loading directions, the pulse is able to effectively scan the ent ire surface of the specimen. 
As soon from the plots the finite element, predictions match well with the theoretical arrival 
times. As tins load and the geometry is symmetric (lie A scan plot at the three source 
locations is the same as seen from Fig(4.7) and Fig(4.8). As seen from the snapshots the 
wave speed in the axial direction is found to be 3755.55 m/sec, which is equal to the shear 
speed of the wave. The wave speed in the circumferential direction is found to be 2545.45 
m/sec. Fig(4.7) and Pig(4.8) contains the A scan plots at the three source locations. 

It is seen from the A-scan plots in the circumferential direction that the wave doesn’t 
attenuate after the completion of the application of the load because it is the surface wave 
which is propogating in the circumferential direction, this is in contrast to the pulse in the 
axial direction which attenuates soon after the completion of the application of the load. 
The reflected pulse from the other end of the tube is clearly seen in Fig(4.7). In Fig(4.8) 
the reflected pulse cannot be observed clearly as the magnitude of the reflected surface 
wave is too small in comparision to the input pulse. 


4.4 Wave interaction with a Circumferential crack 

The objective here is to simulate wave propogation in tubes with circumferential 
cracks and suggest a method to determine the axial and angular location of the crack. The 
simulations are performed on two specimens with different angular and axial locations. 
Fig(4.9) and Fig(4.16) give the geometry of the specimen, loading as well as the location of 
the crack. In order to analyse the plot of and hence to find the location of the crack, A-scan 
plot of differences of displacements at the three locations has been attempted. This be- 
comes necessary as the magnitude of the reflected pulse is much smaller than the magnitude 
of the original pulse, it is observed that the magnitude of the reflected pulse is in the order 
of 10~ 8 , when compared to the original pulse which is of the order of 10 -5 . These plots give 
clear information about the location of the crack. As load is symmetric the A-scan plot of 
from all the source locations should be exactly the same unless there is some flaw which 
creates unsymmetry in the geometry. From the known speed of the propogation of the wave 
and the time taken for the reflection from the crack(flaw), axial location is determined for 
various cases separately. It should be noted that the time obtained from the difference 
plots is the time taken by the wave to travel from the source to the flaw and back to the 
source. From the difference of the A-scan plot one can find angular location of the crack too. 
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4.4.1 Case(i): Crack located at 175 mm length and 180° angle 

Fig(4.9) gives the geometry of the specimen, location and orientation of the crack 
and the loading. The snapshots of the simulation are shown in Fig(4.10) and Fig(4.11). 
Fig(4.12) shows the A-scan plots of displacement difference at all the three source.locations. 
Fig(4.10) and Fig(4.11) shows the propogation of the wave. Fig(4.12) and Fig(4.13) give 
the A-scan plots at all the source locations. Fig(4.14) which give the A-scan plot of dif- 
ferences of(uj — Vj) at various source locations. From these plots it can be seen that the 
displacement suddenly becomes significant at time 137 sec. This corresponds to the time 
at which the reflected wave from the crack surface reaches the source. From the obtained 
time of travel, the axial location of the crack can be obtained. Table(4.2) below gives the 
details of the simulation as given by the axial difference plots . If the A-scan plots (u,- — Vj ) 
difference are closely examined, first plot which shows difference (iq — v 2 ) and the second 
plot which shows the difference (u 2 - v 3 ) are similar, i.e, they are in phase and almost of 
same magnitude. This depicts that angular location of the crack is some where in between 
those of source 1 and 3 along the direction of wave travel(i.e along the direction of loading) 
and preferrably equidistant from both of them, same is the case with plots three and one. 


Table(4.2) 


Time of travel of 

Axial 

Anngular 

Axial 

Angular 


Location 

Location 

Location 

Location 

reflected pulse (/jl sec) 

(estimated) 

(estimated) 

(calculated) 

(calculated) 

137 

175 

165-195 

174.36 

170-190 


4.4.2 Case(ii): Crack located at 233.3 mm length and 225° angle 

Fig(4.16) gives the geometry of the specimen, location and orientation of the crack 
and the loading. The snapshots of the- .simulation are shown in Fig(4.17) and Fig(4.18). 
Fig(4.21) shows the A-scan displacement difference plots at all the three source locations. 
Fig(4.17) and Fig(4.18) shows the propogation of the wave. Fig(4.19) and Fig(4.20) give 
the A-scan plot of displacements at the three source locations. Fig(4.21) gives the difference 
of A-scan plot of(uj - vj) difference at the source locations in the circumferential direction 
from which it can be clearly seen that the displacements becomes significant at time 175 
p sec. This is the time at which the reflected wave from the crack surface reaches the 
source. From the obtained time of travel, the axial location of the crack can be obtained. 
Table (4. 3) below gives the details of the simulation as given by the axial difference plots, to 
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find some information about the angular location A-scan plot of displacement differences 
are to be observed closely. Plots three and one are in the same direction with a small phase 
difference. They have much variation in there peak values, plot three having higher peak 
value than that of plot one. This shows that the angular location of crack is between source 
location 3 and 1 and closer to source location 3. 


Table (4. 3) 


Time of travel of 

Axial 

Anngular 

Axial 

Angular 


Location 

Location 

Location 

Location 

reflected pulse (p sec) 

(estimated) 

(estimated) 

(calculated) 

(calculated) 

175 

233.33 

210-240 

222.72 

230-270 


4.5 Wave interaction with an Axial crack 

The objective here is to simulate wave propogation in tubes with axial cracks and 
determine the axial and angular location of the crack. In case of axial crack the reflections 
are sharp as compered to the Circurn ferential cracks. The simulations are performed on 
two specimens with varying axial and angular locations. 


4.5.1 Case(i): Crack located at 175 mm length and 180° angle 

Fig(4.23) gives the geometry of the specimen, location and orientation of the crack and 
the loading. The results of the simulation are shown in Fig(4.24) and Fig(4.25). Fig(4.26) 
shows the A-scan displacement difference plots at all the three source locations. Fig(4.24) 
and Fig(4.25) shows the propogation of the wave. Fig(4.26) and Fig(4.27) give the A-scan 
plot of the difference at the three source locations. Fig(4.27) give the difference of A-scan 
plot of at the source locations of ( v, — vj) from which it can be clearly seen that the dis- 
placements becomes significant at time 135 p sec. This is the time at which the reflected 
wave from the crack surface reaches the source. From the obtained time of travel, the axial 
location of the crack can be obtained. Table(4.4) below gives the details of the simulation 
as given by the axial difference plots. In order to find some information about the angular 
location A-scan plot of displacement differences need to be observed closely. In Fig(4.27) 
plots two and three are in the same direction with a small phase difference. This depicts 
that the crack is between source 2 and 1 and close to source location 2. It is also observed 
that the plot one is exactly 180 degrees apart from that of two and three, this means that 
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the crack is very close to source two. 


Table(4.4) 


Time of travel of 

reflected puls e(/i sec) 

Axial 

Location 

(estimated) 

Anngular 

Location 

(estimated) 

Axial 

Location 

(calculated) 

Angular 

Location 

(calculated) 

135 

175 

180 

171.81 j 

170-190 


4.5.2 Case(ii): Crack located at 140 mm length and 120° angle 

Fig(4.30) gives the geometry of the specimen, location and orientation of the crack 
and the loading. The snapshots of the simulation are shown in Figs(4.31) and (4.32). 
Fig(4.35) shows the A-scan displacement difference plots at all the three source locations. 
Fig(4.33) and Fig(4.34) shows the propogation of the wave. Fig(4.33) and Fig(4.34) give 
the A-scan plot of at the three source locations. Fig(4.35) gives the difference of A-scan 
plot of (v{ — Vj) source locations in the circumferential direction from which it can be clearly 
seen that the displacements suddenly became significant at time 135 /.i sec. This is the time 
at which the reflected wave from the crack surface reaches the source. From the obtained 
time of travel, the axial location of the crack can be obtained. Table(4.5) below gives the 
details of the simulation as given by the circumferential difference plots. In order to find 
some information about the angular location A-scan plot of displacement difference are 
to be observed closely. Plots three and one are in the same direction with a small phase 
difference. This may depict that the crack is between source 3 and 2 and close to 2. 


Table(4.5) 


Time of travel of 

Axial 

Anngular 

Axial 

Angular 


. Location 

Location 

Location 

Location 

reflected pulsed sec) 

(estimated) 

(estimated) 

(calculated) 

(calculated) 

I— 1 

O 

o 

140 

120 

128 

100-140 
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FIG 4.9: SPECIMEN GEOMETRY AND LOADING 
CIRCUMFERENTIAL CRACK: L=I75mm ,0-?iSO deg 
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FIG 4.1 l(cont.): SNAPSHOTS OF AXIAL DISPLACEMENTS 

CIRCUMFERENTIAL CRACK: L=175mm ,0-180 deg 
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FIG 4.15: A-SCAN PLOTS OF DIFFERENCE OF CIRCUMFERENTIAL DISPLACEMENT 

AT THE SOURCE LOCATIONS 
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FIG 4.22: A-SCAN PLOTS OF DIFFERENCE OF CIRCUMFERENTIAL DISPLACEMENT 

AT THE SOURCE LOCATIONS 



FIG 4.23: SPECIMEN GEOMETRY AND LOADING 
AXIAL CRACK: 1. 1 75mm, ISO tit'!* 
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FIG 4.24(cont.): SNAPSHOTS OF CIRCUMFERENTIAL DISPLACEMENTS 
AXIAL CRACK: L= 175mm, 0=180 deg 
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FIG 4.30: SPECIMEN GEOMETRY AND LOADING 
AXIAL TRACK I . I lOium, i~> I ’it dot* 
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FIG 4.31(cont.): SNAPSHOTS OF CIRCUMFERENTIAL DISPLACEMENTS 
AXIAL CRACK: L= 140mm, 0=120 deg 
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Chapter 5 
Conclusions 


5.1 Conclusions 

A finite element model has been developed to study the elastic wave propogation 
and scattering in thin tubes. A number of simulations are performed to suggest a method 
to determine the axial and angular location of the cracks in thin tubes. The finite element 
method has been found to be very effective in describing the complex wave/defect interac- 
tion. 

From the numerical computation point of view it is concluded that the Newmark’s time 
integration scheme along with degenerated shell elements is optimum combination to solve 
the hyperbolic partial differential equation governing the wave motion in tubes. 

From the results of the simulations performed the following conclusions are drawn: 

1. A methodology has been suggested to for the detection of location and classification 
of defects in thin tubes. 

2. The method developed helps in clearly locating the axial location of the defect. 

3. This method qualitatively gives the angular location of the defect. 

4. Using this method long length of the tube can be inspected from each transducer. 
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5.2 iScopo of fufnro work 


LDik! to computational ms I, notions only a small domain was modelled in the present work. 
I'uthei analysis can be carried out with a liner mesli to get more realistic results. 


2. In this work angular location of the crack is found qualitatively, more realistic results 
can be found by making a quantitative analysis. 

3. Transducer models can be incorporated into the code so that the results can be closer to 
the practical inspection condition. 
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